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Motivated by the idea that some characteristics are specific to the relations between individuals 
and not of the individuals themselves, we study a prototype model for the dynamics of the states 
of the links in a fixed network of interacting units. Each link in the network can be in one of two 
equivalent states. A majority link-dynamics rule is implemented, so that in each dynamical step 
the state of a randomly chosen link is updated to the state of the majority of neighboring links. 
Nodes can be characterized by a link heterogeneity index, giving a measure of the likelihood of 
a node to have a link in one of the two states. We consider this link-dynamics model on fully 
connected networks, square lattices and Erdos-Renyi random networks. In each case we find and 
characterize a number of nontrivial asymptotic configurations, as well as some of the mechanisms 
leading to them and the time evolution of the link heterogeneity index distribution. For a fully 
connected network and random networks there is a broad distribution of possible asymptotic 
configurations. Most asymptotic configurations that result from link-dynamics have no counterpart 
under traditional node dynamics in the same topologies. 
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I. INTRODUCTION 

Collective properties of interacting units have been 
traditionally studied considering that each of these units 
has a property or state, and interacts with others in a 
network of interactions. The result of the interaction 
depends on the state of the interacting units. For 
example, in a spin system in a lattice, there is a spin 
in each node with a given state and it interacts with 
its neighbors in the lattice, in a way that depends on 
their relative spin state. The same basic set up has 
been implemented in individual or agent based models 
of social collective properties [l|. These models endow 
individuals with a variable, which can be discrete or 
continuous, describing for example, an opinion state. 
The model also prescribes a dynamical rule, which 
results in changes of the states of the agents that depend 
on the state of the agents with whom they interact. 
However, there are a number of possible interactions 
among individuals in which the state variable is more 
properly described as a state of the interaction link 
than a state of the individuals in interaction. This is 
specially the case for relational interactions such as 
friendship, trust, method of communication (phone or 
skype), method of salutation (kiss or handshake), etc. 
It is also the case in language competition dynamics 
So far language has been modeled in this context 
as an individual property, but use of a language, 
as opposed to knowledge of a language, is more a link 



than a node property in a social network of interactions. 
Noteworthy, data on link states associated with trust, 
friendship or enmity, obtained from on-line games and 
on-line communities, is now available and has been 
recently analyzed [f| p . 

Social balance theory Q is a well established prece- 
dent in the study on link states and link interactions. 
This theory states that a relation can be positive or 
negative and that there is a natural tendency to avoid 
unbalanced triads of individuals, where an unbalanced 
triad is defined as one for which the multiplication of the 
states of the three links is negative. A number of recent 
studies address social balance in complex networks, im- 
plementing stochastic link dynamics that explore when 
a balance situation is or it is not reached asymptotically 
Social balance theory has also been confronted 
with large scale data 0, Hf, and alternative theories 
for the interaction of positive/negative relations have 
also been proposed @, . Focusing on link properties 
has also been emphasized in the problem of community 
detection in complex networks |12rfla | . This opposes the 
traditional view of identifying network communities with 
a set of nodes [l7j . and it makes possible for an individ- 
ual to be assigned to different communities. Finally, the 
idea of considering link dynamics is also present in the 
problem of network dynamics controllability 18[. Here 
the aim is to identify the most relevant links to drive the 
system to a desired global state of the network, instead 
of focusing on the dynamically most influent nodes (l9j | . 
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example in social balance, where the two link states play 
different roles. Equivalent link states can occur in many 
relational interactions including, for example, salutation 
or competition of languages of the same prestige. As 
a first step towards the characterization of such link 
dynamics we choose to investigate a majority rule 
dynamics akin to a zero-temperature kinetic Ising model 
but for the states of the links, instead of the state of the 
nodes. We show that such link majority rule dynamics 
on complex networks results in a variety of non-trivial 
different asymptotic states which are generally not found 
when studying traditional node-dynamics in the same 
topologies. We also show how a quantity characterizing 
the node behavior naturally arises for the link states, so 
that nodes can also be characterized by the state of the 
links connected to the node. 

The paper is organized as follows: in section II we 
define our majority rule link dynamics model, as well as 
some quantities introduced for its characterization. In 
sections III, IV and V we describe our results on a fully 
connected network, a square lattice and Erdos-Rcnyi 
random networks, respectively. Section VI contains a 
discussion summary. Besides the material presented in 
this paper, which is self-explanatory, the supplementary 
material available in the web will help visualize different 
aspects of the model on different networks. 



II. A MODEL FOR LINK DYNAMICS 



its links dynamics 
kii Degree of node i. 

lf^ B - Number of A/B edges connecting node i. 

p: Density of nodal interfaces: It serves as order param- 
eter. It gives a local measure of the level of order 
in the system. It vanishes when the system is com- 
pletely ordered. 

V^JV ]A]B 
_ Z-ii=l H H 

~E£ifc(*i-i)/2' 

be. Link heterogeneity index of node i. It is a node char- 
acterization giving a measure of the likelihood of a 
node to have A or B links (6, = +1 or bi = —1 for 
all links of the same type, bi = for a completely 
symmetric case). 



P(b,t): Link heterogeneity index distribution, probabil- 
ity that a randomly chosen node has link hetero- 
geneity index b at time t. 

S(t): Survival probability, probability that a realization 
of the stochastic link dynamics has not reached a 
fully ordered configuration at time t. 



We consider a fixed network composed by N nodes and 
L edges. The state of each link i — j is characterized by a 
binary variable Si-j which can take two equivalent values 
A ov B. We study a majority rule for the dynamics of 
the state of the links. The basic step of the dynamics is 

i. Randomly choose a link i — j. 

ii. Update its state to the one of the majority of links 
in its first neighborhood (two links are considered 
first neighbors if they are attached to a common 
node). In case of a tie, the state of the link is 
randomly chosen 

The time unit is set to TV basic steps so that for each node, 
on the average, the state of one of its links is updated per 
unit time. 

There exist two trivial absorbing ordered configura- 
tions, for which all the links in the system have the 
same state. The dynamics tend to order the system 
locally. We investigate whether, depending on the 
topology of the network, the dynamics orders the system 
globally or if the system reaches asymptotic disordered 
configurations with coexistence of both states. We will 
also characterize the transient dynamics towards these 
asymptotic states. For these purposes we will consider 
the following quantities characterizing the network and 



III. FULLY CONNECTED NETWORK 

We first consider a fully connected network in which 
every node is connected to every other node so that L = 
N(N — l)/2. Note however that every link is not a first 
neighbor of every other link. 

A. Time evolution 

Fig. Q] shows the time evolution the ensemble average 
of the order parameter < p > and the survival proba- 
bility S(t) for random initial conditions. The average 
order parameter shows a decay towards a plateau, 
indicating that the absorbing ordered configurations are 
not always reached. Comparing this result with the 
survival probability, which also saturates at a certain 
value after a transient, we conclude that the plateau in 
the average order parameter is due to realizations which 
get frozen in a configuration with coexistence of states. 
The analysis of single realizations of the link dynamics 
(lower panel of Fig. [1} shows smooth dynamics to an 
asymptotic state in which the order parameter is frozen. 
In the following we investigate the characteristics of 
these frozen asymptotic configurations. 
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FIG. 1. Upper panel: Evolution of the average order parame- 
ter on a fully connected network. Inset: Survival probability. 
TV = 100 for the black solid line, TV = 300 for the red dashed 
line and TV = 600 for the blue dashed-dotted line. Averages 
are taken over 10 3 realizations. Lower panel: Evolution of 
the order parameter for single realizations of the dynamics 
on a fully connected network of size TV = 300. We show two 
different kinds of realizations: a realization reaching an ab- 
sorbing ordered state (solid line) and a realization ending in 
a disordered frozen configuration (dashed line). 



B. Asymptotic states 

The probability of having a certain value of p in the 
asymptotic configurations is plotted in Fig. [2] We ob- 
serve a very heterogeneous set of possible final configura- 
tions in addition of the most probable ordered configura- 
tion (p = 0). The disordered frozen configurations can be 
classified by the number n& of different link heterogene- 
ity indices occurring in each configuration. The limiting 
case rij, = 1 corresponds to the ordered configuration in 
which all nodes have the same heterogeneity index 6=1 
or b = —1. 



1. Simplest frozen configurations 

The simplest frozen configurations on a fully con- 
nected network are of the type shown in Fig. Oa). 
They consist of a set of k nodes that have only links 
in one state (red links), and the rest of nodes, TV — k, 
having all their links in the other state (blue links), 
except for the links with the k nodes of the first set. 
It is clear that for this kind of configuration there are 
only two types of nodes in terms of link heterogeneity 
index. The group of size k having \b\ = 1 and therefore 
contributing to the asymptotic p(b, t = oo) in the peaks 
b = ±1 (see Fig[6] a), and another group of size TV — k 
with |b| = (2k - TV - 1)/(TV - 1). Therefore, for these 
configurations n& = 2. 
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FIG. 2. Probability of having a certain value of the order pa- 
rameter in the asymptotic configuration for a complete graph. 
The calculation is done over 10 4 realizations for system sizes 
TV = 100 (black circles), TV = 300 (red squares) and TV = 600 
(blue diamonds). 
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FIG. 3. a) Simple frozen configurations in a fully connected 
network, b) Frozen configuration with ni = 3 on a fully 
connected network. 



These configurations can be proven to be frozen for a 
range of sizes k. For this purpose one has to find how 
many of the neighboring links of a given link are in each 
state and impose that links in state A (B) have more 
A (B) neighbors than B (A) neighbors. In this way one 
easily concludes that configurations as the one in Fig. [3] a) 
are frozen whenever 

1< k < TV/2 - 1. 

These solutions exist for TV > 4. In Fig. |4] we show the 
probability density to reach a configuration with a cer- 
tain fraction fc/TV with |6| = 1, given that the asymptotic 
configuration is of the type shown in Fig. [3]a). All the 
possible configurations can be reached from random ini- 
tial conditions. 



4 



10 1 




10 0.1 0.2 0.3 0.4 0.5 
k/N 

FIG. 4. Probability density of getting to a simple frozen con- 
figuration like the one in Fig. 0a) with a certain fraction k/N 
of nodes with \b\ = 1, starting from random initial condi- 
tions on a complete graph. Sizes are N = 100 (black circles), 
N = 300 (red squares) and N = 600 (blue diamonds). The 
statistics are over 10 realizations of the system. 




FIG. 5. Probability of reaching a frozen configuration with 
a certain number of different link heterogeneity indices rib, 
starting from random initial conditions on a complete graph. 
Sizes are N = 100 (black circles), N = 300 (red squares) 
and -/V = 600 (blue diamonds), and the statistics are over 10 5 
realizations of the system. 



2. Other asymptotic configurations 

Figure [5] shows the probability of reaching an asymp- 
totic configuration with a certain number of different 
link heterogeneity indices n& m the system. The or- 
dered configurations rib = 1 and the ones with rib = 2 
described above are the most probable. An example 
of a configuration with rib = 3 is shown in Fig. [3jb). 
These configurations have k nodes with |6| = 1, I nodes 
with b = (2k — N + 1)/(N - 1) and N - k - I nodes 
with b = (N - 21 - \)/{N - 1). Following the same 
argument used for rib — 2 configurations, we can con- 
clude that rib — 3 configurations are frozen provided that 



k > 1 

I < N/2 - 1 
k < N/2 - 1 
I > k + 1 

When rib is increased, the frozen configurations become 
structurally more complicated, and are much less 
probable. Empirically we have found that always a 
group of agents with |6| = 1 appears. To characterize a 
frozen solution with rib we need rib ~ 1 parameters and 
rib{rib + l)/2 inequalities, which arise imposing that the 
state of every link is frozen and give a boundary for the 
possible architecture of those configurations. 



C. Link heterogeneity index distribution 

Figure [51a shows the time evolution of the link het- 
erogeneity index distribution: We observe that it evolves 
from a distribution peaked around 6 = for random ini- 
tial conditions, to a bimodal distribution peaked around 
b = ±1, with a quite homogeneous probability of hav- 
ing any link heterogeneity index. This statistics char- 
acterized by this distribution includes the most proba- 
ble realizations that reach ordered states but also other 
which freeze in configurations with nodes with different 
possible value of b, as discussed in the characterization 
of the asymptotic configurations. Note that both type 
of realizations contribute to the peaks at b = ±1 since 
frozen disordered asymptotic configuration have at least 
one group of agent with b = ± 1 . 



IV. SQUARE LATTICE 

In order to account for local interaction effects we first 
consider a square lattice with periodic boundary condi- 
tions. 



A. Time evolution 

Fig. [7] shows the time evolution the ensemble average 
of the order parameter < p > and the survival prob- 
ability S(t) for random initial conditions. We observe 
a qualitative behavior very similar to the one in Fig. Q] 
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FIG. 6. Probability density distribution of link heterogeneity 
index P(b,t) for different times averaged over 10 3 realizations. 
The initial condition is in black circles. Time evolution is in 
the following order: red squares, green diamonds, blue up tri- 
angles and magenta left triangles, a) Fully connected network 
of size N = 100 at times 0, 50, 100, 200, 500. b) Square lattice 
of size N = 2500 for times 0, 500, 1000, 2000, 3000. In this case 
the actual probability of having index b is plotted and not the 
density, c) Erdos-Renyi random network of size TV = 1000 and 
average degree (k) = 10 at times 0, 50, 100, 200, 500. All plots 
are symmetric around 6 = due to the equivalent nature of 
the states A and B. 
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FIG. 7. Upper panel: Evolution of the average order parame- 
ter on a square lattice. Inset: Survival probability. N = 2500 
for the black solid line, N = 3600 for the red dashed line 
and N = 4900 for the blue dashed-dotted line. Averages are 
taken over 10 3 realizations. Lower panel: Evolution of the 
order parameter for single realizations of the dynamics on a 
square lattice of size TV = 2500. We show three different real- 
izations, corresponding to the three possible asymptotic con- 
figurations: ordered state (dashed line), vertical/horizontal 
single stripe (solid line) and diagonal single stripe (dotted- 
dashed line). 



for a fully connected network, i.e. (p) and S(t) decay 
smoothly to a plateau value. Together with the plot 
of single realizations of the stochastic dynamics (lower 
panel of Fig. [7] ) this indicates that some of the realiza- 
tions reach an asymptotic ordered state, while others get 
trapped in a disordered configuration for which the order 
parameter remains constant for all times. We have found 
only three different types of realizations characterized by 
their asymptotic configurations, as we discuss next. 



B. Asymptotic states 

The probability of reaching one of the three main pos- 
sible asymptotic states, characterized by their value of 
the order parameter, is shown in Fig. [5] These configu- 
rations are depicted in Fig. [9] Note that these are the 
asymptotic configurations starting from random initial 
conditions and using periodic boundary conditions. 



• Ordered configurations: 
state and p = 0. 



All links are in the same 



• Dynamically trapped configurations, where the or- 
der parameter remains constant, p = 1/y/N, but 
the densities of links in each state fluctuate around 
a certain value. These configurations form verti- 
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FIG. 8. Probability of reaching a given asymptotic value of 
the order parameter on a square lattice with periodic bound- 
ary conditions and starting from random initial conditions. 
There are three different possible configurations, namely or- 
dered state, horizontal/vertical stripes and diagonal stripes. 
Sizes are N = 2500 (black circles), N = 3600 (red squares) 
and N — 4900 (blue diamonds). Statistics computed from 
10 4 realizations. 
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FIG. 9. Different asymptotic disordered configurations on a 
square lattice with periodic boundary conditions, a) Verti- 
cal/horizontal single stripe. The gray links keep changing 
state forever, while all other links are in a frozen state, b) 
Diagonal single stripe. All links are frozen, c) Percolating 
diamond. All links are frozen. 



cal/horizontal stripes, as shown in Fig. [9^,. These 
configurations are dynamical traps from which the 
system cannot reach the ordered state: links in the 
boundaries of the stripes continue to blink without 
changing the value of the order parameter. Sin- 
gle stripe configurations are the ones reached from 
random initial conditions, but configurations with 
a larger number of stripes and a value of the order 
parameter which is a multiple of l/y/N are also 
dynamical traps of the model. 

• Frozen configurations, where the order parameter 
and the densities of links in each state remain con- 
stant. Configurations reached from random ini- 
tial conditions are single diagonal stripe as the one 
shown in Fig. [!Tp, with a value of the order pa- 
rameter p = 3 1^ ■ There are other frozen con- 
figurations for our dynamical model which we have 
not observed in our simulations with random initial 
conditions. These include multiple diagonal stripes 
and a combination of diagonal front that we call 
percolating diamond (Fig.[5p): It contains a square 
of links in one state, rotated an angle of 45 degrees, 
surrounded by links in the opposite state and which 
percolates through the network. For the percolat- 
ing diamond configuration p = ^ . 



C. Link heterogeneity index distribution 
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FIG. 10. Upper panel: Evolution of the average order pa- 
rameter on Erdos-Renyi networks of average degree (k) = 10. 
N = 1000 for the black solid line, N = 5000 for the red 
dashed line and N — 10000 for the blue dashed-dotted line. 
Averages are taken over 10 3 realizations of different initial 
conditions and different realizations of the random network. 
Lower panel: Evolution of the order parameter for single re- 
alizations of stochastic dynamics on an Erdos-Renyi random 
network of size N — 1000 and average degree (k) = 10. Two 
different realizations are shown, each one ending in a different 
configuration with frozen order parameter. 



vertical/horizontal single stripe. For these configurations 
there are nodes whose heterogeneity index keeps 

jumping from b = ±1 and b = ±0.5. Last, the probability 
of having nodes with 6 = comes from the third possible 
asymptotic configuration, diagonal single stripe. In this 
configuration 2y/~N — 2 nodes have 6 = and the other 
nodes are divided into two equal groups with 6 = 1 and 
6= -1. 



V. RANDOM NETWORKS 

In order to account for the role of network heterogene- 
ity we finally consider the link dynamics model on Erdos- 
Renyi random networks. 



For a square lattice the link heterogeneity index takes 
values 6 = ±1, ±0.5, 0. The evolution of the distribution 
P(b,t) is shown in Fig. ®b). It evolves from an initial 
peak at 6 = to a distribution with two peaks at 6 = ±1 
and a minimum value at 6 = 0. This evolution can be 
understood from the asymptotic configurations described 
above: The two peaks at 6 = ±1 originate in the most 
probable ordered configurations, but also on the large 
percentage of nodes with 6 = ±1 in the two other possible 
asymptotic states. The values 6 = ±0.5 appear only 
in the second most probable asymptotic configuration: 



A. Time evolution 

Proceeding as in the previous cases we show in Fig. [10] 
the time evolution of the ensemble average order param- 
eter. The survival probability (not shown) is one at all 
times, except for small systems or networks of high aver- 
age degree where it tends to a fully connected like behav- 
ior. Our results indicate that all stochastic realizations 
of the dynamics reach an asymptotic disordered configu- 
ration with a constant value of p. 
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FIG. 11. Probability of having a certain value of the or- 
der parameter in the asymptotic configuration on a random 
graph. The calculation is done over 10 4 realizations for system 
size N = 1000 and average degrees (fe) = 10 (black circles), 
(k) = 20 (red squares) and (k) = 40 (blue diamonds). 



Frozen links 

FIG. 12. Example of change in state which changes the den- 
sities of blue and red links conserve the value of the order 
parameter p. Independently of the state of the grey link this 
motif will contribute to the order parameter of the whole sys- 
tem with p — 1/5. 



links, i.e. links which do not change state. Ini- 
tial conditions and the particular topology of the 
realization will determine how frequent is this ef- 
fect and whether this leads or not to an ordered 
configuration. 



B. Asymptotic states 

We observe a large variety of asymptotic states char- 
acterized by different values of the order parameter p, as 
shown in Fig. QTJ Increasing the average degree of the 
network, the distribution of final values of p approaches 
the one for a fully connected network (Fig. [2]): The dis- 
tribution develops a peak that moves towards p = and 
another peak near the maximum possible asymptotic or- 
der parameter value. 

For random networks we also find three kinds of 
asymptotic states: 

• Ordered configurations: All links are in the same 
state and p = 0. This configuration is only observed 
in small systems or in systems with high average 
degree (close to fully connected network). 

• Dynamically trapped configurations: the order pa- 
rameter remains constant, but the densities of links 
in each state vary with time. 

• Frozen configurations, where the order parameter 
and the densities of links in each state remain con- 
stant. 

It is possible to identify some basic mechanisms leading 
to the observed traps. Among them: 

• The role of hubs: If a node i is such that hi 3> kj 
for any neighboring node j, then the links attached 
to node i usually end up sharing all the same state, 
which in most cases is the one of the initial majority 
state in that set of links. This effect creates frozen 



• Dynamics conserving the value of p: There exist 
changes of the state of the links which do not cause 
a change in the value of p. These changes are those 
for which the link changing state has a symmetric 
environment, with the same number of neighbors 
in each state as shown in Fig. Q21 This situation is 
the one also found in a square lattice Fig. [!|]a. This 
kind of phenomenon can appear in more complex 
forms, as shown in Fig. 1131 There one can see that 
the order parameter is frozen after approximately 
10 time steps, but the configuration of the system 
keeps changing, as can be seen from the snapshots 
of the system configuration at different times. 



C. Link heterogeneity index distribution 

The evolution of the distribution of link heterogene- 
ity indices in random networks is shown in Fig. [S]c). 
The initial distribution is broad, but smoothly peaked 
around 6 = 0. This evolves to a bimodal distribution 
peaked around b = ±1. The hubs of the network are 
prone to become nodes with b = ±1, which in turn pulls 
more nodes to this value of b. The fact that links can be 
in frozen states for different parts of the network implies 
that between patches of ordered domains, there are nodes 
with any value of the link heterogeneity index. This con- 
tributes to the broad distribution of b values between the 
two peaks. Blinking links in dynamically trapped con- 
figurations also contribute to the broad distribution of 
intermediate values of b. 
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FIG. 13. One realization on a small random network of size 
N = 20. Top left panel shows the evolution of the order pa- 
rameter, which freezes after approximately 10 time steps. The 
other panels show the configuration of the system at different 
times. The color of the nodes reflects their link heterogeneity 
index. Red (blue) is for having all links in the red (blue) op- 
tion, white is for having half of the links in each color. The 
changes in the configuration do not affect the value of the or- 
der parameter. For example the only difference between the 
configuration at t = 20 and the one at t = 120 is the state 
of a single link. If we count we can see that the link has the 
same number of neighbors in each state. One can check that 
all the changes of state are of the type depicted in Fig. [12] 



latter do not have an analog in the corresponding node 
dynamics. In the case of Erdos-Renyi random networks 
we have described the mechanisms leading to the exis- 
tence of very heterogeneous asymptotic configurations 
which are either frozen or dynamical traps. 

It is clear that a link-dynamics model can be mapped 
into an equivalent node-based problem by changing the 
network of interaction. The node-equivalent network 
is the line-graph [2fjl - [22j of the original network. The 
line-graph is a network where the links of the original 
network are represented by a node and are connected to 
those nodes that represent links that were first neighbors 
in the original network. This mapping of the problem 
has not been pursued here since it obscures our original 
motivation and, given the complexity of the line graphs 
of the networks considered here, it has been found not 
to be particularly useful for a quantitative description 
of the dynamics. However, the mapping does provide 
additional qualitative understanding of our findings: 
The line graph is a network with higher connectivity 
since all links that converged originally in a node form 
a fully connected subgraph in the line-graph, as clearly 
seen in the line-graph of a fully connected network or 
a square lattice. This results in an increased cliqueness 
of the line graph, as compared to the original network. 
Such cliqueness is behind the topological traps that give 
rise to the wide range of possible asymptotic states that 
we find for the link-dynamics. In addition, the mapping 
of a hub of the original network in the corresponding 
line-graph also helps understanding the different role 
played by hubs in node or link-based dynamics: as 
discussed in Section IV, hubs tend to freeze link states 
in their neighborhood. 



VI. SUMMARY AND DISCUSSION 

The study of a majority rule for the dynamics of 
two equivalent link states in a fixed network uncovers a 
set of non-trivial asymptotic states which are generally 
not present when studying the classical node-based 
majority rule dynamics. The characterization given of 
the asymptotic states in fully connected networks, square 
lattices and Erdos-Renyi random networks provides the 
basis for the understanding of the evolution of the link 
heterogeneity index distribution. For a fully connected 
network and for a square lattice we have fully character- 
ized the asymptotic states reached from random initial 
conditions. In a fully connected network we have found 
large heterogeneity in the asymptotic configurations. 
All these configurations, classified by the number rib of 
heterogeneity indexes present in the configuration, are 
frozen. Note that for the corresponding node-dynamics 
in the same network only an asymptotic ordered con- 
figuration is found (rib — !)• In a square lattice we 
have found asymptotic configurations which are ordered, 
frozen and disordered, or dynamically trapped. The 



The link heterogeneity index is a useful way of charac- 
terizing nodes in a given link configuration. For example 
in node based models of language competition, a node 
can be in state A or B corresponding to two competing 
languages, and bilingualism can only be introduced 
through a third node bilingual state AB In the 

framework of link dynamics, state A or B characterizes 
the language used in a given interaction between two 
individuals, and the link heterogeneity index is a natural 
measure of the degree of bilingualism of each individual 
(node). Continuing with this example, a next step is to 
consider the mixed dynamics of language competence 
(node dynamics) and language use (link dynamics). In 
general, consideration of the coevolution of link and 
node states is a natural framework that emerges in the 
study of collective behavior of interacting units. In 
physical terms, the states of the interacting particles 
are coupled to the state of the field that carries the 
interaction. Another possible avenue of research is the 
addition of more realistic features to the model, such as 
the temporal patterns of human interactions (2a , |2~H . 
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